Quasiparticle and Cooper Pair Tunneling in the Vortex State of Bi 2 Sr 2 CaCu 2 08+5 
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From measurements of the c-axis FV characteristics of intrinsic Josephson junctions in 
Bi2Sr2CaCu208+i (Bi-2212) mesas we obtain the field dependence (H || c) of the quasiparticle 
(QP) conductivity, a q (H,T), and of the Josephson critical current density, J C (H,T). o q {H) in- 
creases sharply with H and reaches a plateau at 0.05 T< H < 0.3 T. We explain such behavior by 
the dual effect of supercurrents around vortices. First, they enhance the QP DOS, leading to an 
increase of a q with H at low H and, second, they enhance the scattering rate for specular tunneling 
as pancakes become disordered along the c-axis at higher H, leading to a plateau at moderate H . 



PACS numbers: 74.25.Fy, 74.50. +r, 74.72.Hs 

Important information regarding pairing symmetry in 
cuprate superconductors has been obtained from the 
quasiparticle (QP) spectrum. Strong evidence in favor of 
d-wave symmetry of the superconducting order parame- 
ter has come from ARPES and STM experiments |],|. 
The effect of a magnetic field on QP properties becomes 
particularly interesting because, in the nodal regions of 
a d-wave gap, it probes fine details of the low-energy QP 
spectrum and its alteration by supercurrents surround- 
ing vortices. An increase of QP density of states (DOS) 
produced by the Doppler shift of near-nodal QP ener- 
gies by the vortex supercurrents is predicted to lead to 
an increasing QP DOS with magnetic field HQ. This 
effect was apparently observed in calorimetric measure- 
ments in HTS p). In-plane thermal conductivity mea- 
surements however have indicated a more complex mag- 
netic field dependence. At low temperatures K a b(H) was 
seen to increase sublinearly with field in agreement with 
predictions of the Volovik mechanism [pi . At higher tem- 
peratures, T >2-5 K. K a b{H) was observed in Bi-2212 
to first decrease with increasing magnetic field and then 
reach a plateau and remain constant up to fields of several 
tesla ||. A similar behavior was observed also in YBCO 
. Franz || attempted to explain the high temperature 
plateau behavior as resulting from a compensating effect 
of QP scattering from in-plane vortex disorder. While 
this provides a qualitative description of the high tem- 
perature behavior, it fails to account for the fact that 
the low temperature sublinear increase with H shows no 
effect of a compensating scattering from vortex disor- 
der. Thus, the thermal measurements have presented an 
ambiguous picture of the competing contributions of vor- 
tices, increasing DOS and scattering, to QP transport. 

More direct information on QP properties can be ob- 
tained by study of the charge transport, namely, the 
c-axis conductivity a c (H), in Josephson-coupled lay- 
ered HTS. Recent measurements on micron-sized mesas 
of Bi-2212 have demonstrated that QP conductivity in 
the superconducting state can be accessed by exceeding 
the Josephson critical current, switching the interlayer 



Josephson junctions into the resistive state and return- 
ing to low currents on the resistive branch |Tc|| . Recent 
theoretical development has shown that c-axis QP con- 
ductivity is determined by the QP DOS and the effective 
scattering rate for interlayer tunneling [pd] . In this Letter 
we provide an unambiguous determination of the effects 
of vortices on the QP DOS and on the effective scattering 
rate for c-axis tunneling. 

We study the effect of vortices on the c-axis QP con- 
ductivity, tj q (H,T), by measuring the I-V characteristics 
in the resistive state of small Bi-2212 mesas as a func- 
tion of the magnetic field < H < 9 T applied parallel 
to the c-axis. At the same fields we extract the depen- 
dence of the Josephson critical current density J C (H,T) 
to obtain independent information on c-axis correlation 
of pancake vortices. Our main finding is that <r q {H,T) 
first increases rapidly with field, then reaches a plateau 
at 40 mT < B < 0.5 T due to the dual role of vortices. 
As the field further increases, cr q (H) increases linearly, as 
was found previously p2] |. We observe a moderate drop 
of J C (H, T) with H in the field range of the plateau and 
then a faster drop of J c in fields where the linear growth 
of <r q (H) takes place. This drop indicates strong disorder 
in vortex positions along the c-axis in high fields in agree- 
ment with enhancement of the effective scattering rate at 
interlayer tunneling responsible for plateau formation. 

For our experiments we used several step-like mesa 
samples with T c — 73 K and area 5 = 2x2 (im 2 manu- 
factured from high-quality Bi-2212 whiskers utilizing the 

TABLE I. Bi-2212 mesa samples; S w 4 ^m 2 , T c w 73 K. 
R n is the resistance at T = 300 K, N is the number of intrinsic 
junctions, J c is the maximum critical current density (H = 0) , 
7 is the in-plane effective scattering rate, and a q (0,0) is the 
QP conductivity at H = 0, T ->• 0. 
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ml 


2.1 


42 


830 


32.0 


1.5 


m2 


4.5 


90 


500 


31.8 


2.1 


m3 


0.85 


17 


1100 


34.0 


4.3 


m4 


3.0 


60 


450 




1.8 
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Focused Ion Beam technique 10 . In this Letter we 
present data obtained on 4 typical samples (See Table 
I). We performed measurements of I-V characteristics at 
different magnetic fields and temperatures in a standard 
He-flow cryostat, which provided temperature stabiliza- 
tion better than ±10 mK. Magnetic fields up to 9 T 
were applied along the c-axis of the mesa by a super- 
conducting solenoid. The intrinsic Josephson junctions 
in high quality mesas have a small shunting conductance 
due to QP tunneling and thus behave as underdamped 
junctions with highly hysteretic I-V dependence. Typi- 
cal multibranching I-V curves in zero-field-cooled (ZFC) 
mode at T = 4 K are shown in Fig. 1. The increasing 
current branch provides information on the critical cur- 
rent, whereas the decreasing current branch corresponds 
to the resistive state, where current is due to QP tunnel- 
ing. Direct comparison of I-V curves at H =0 T, 0.05 T, 
and 0.15 T (Fig. 1) shows that a) critical currents drop 
with H, and thus all branches collapse to that of the 
resistive state, and b) variation of critical current for dif- 
ferent junctions increases with H in comparison with that 
at H = fl 
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FIG. 1. I-V curves for Bi-2212 mesa m3 at H = 0, 
H = 0.05 T, and H = 0.15 T. For this sample all 17 branches 
are resolved. 

It is important that the branch of I-V curves for de- 
creasing current, when all junctions are in the resistive 
state, is the same for all current sweeps for both field cool- 
ing (FC), and ZFC (Fig. 2). This observation allows us to 
study the field dependence of QP conductivity, a q (H, T) , 
measured in ZFC. The value of average QP resistivity, 
(pq), was obtained from the all-junction resistive-state 
I-V curve at low currents as described in Ref. [|l0|, and 
a q = l/(p q ) was calculated. 

In contrast to the behavior of a q (H), the distribution 
of critical current over junctions is different for FC and 
ZFC modes. The distribution of J C (B,T) over junctions 
in small-area mesas arises because a) at T — > and H = 
junctions are inequivalent with respect to J c at least due 
to their geometrical position in the mesa, b) at nonzero 
field J c (B,0) for a junction between layers n and n+ 1 



depends on vortex positions in the layers, which vary 
randomly with n in the case of uncorrelated pinning [ p"3| , 
and c) at nonzero T, jumps from the superconducting 
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FIG. 2. I-V curves for H = 0.03 T, FC (o) and ZFC (▼) 
modes. In the FC mode vortices initially are better correlated, 
providing larger J c . Switching the junctions in the resistive 
state suppresses c-axis pancake correlations, and the system 
becomes similar for both FC and ZFC. Inset: the distribution 
of J c in junctions for ZFC (dash) and FC (line) modes. 

state to a resistive one occur at different J due to thermal 
fluctuations |Q. We can find the distribution function 
of the critical current density, /(J), from our data as 
/(J) = d(Vt/Vi)/dJ, where V T (J) [V](J)] are voltages 
on increasing [decreasing] currents. The inset in Fig 2 
shows smoothed distribution functions f(J) for FC and 
ZFC modes. We note that average critical current den- 
sity, (J C (H = 0.03T)), obtained in the FC mode at 4 K 
is larger by the factor rs 2 than that in the ZFC mode. 
This result is in agreement with measurements fl5f| of the 
Josephson plasma resonance ( JPR) frequency. Below the 
irreversibility line (in the vortex glass phase) the JPR 
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FIG. 3. For mesa ml a q (V — > 0) is plotted vs. H for three 
different temperatures. 
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frequency in the FC mode was found to be larger than 
in the ZFC mode. A very important new observation is 
that curves obtained in the second and subsequent cur- 
rent sweeps in the FC mode are similar, and they prac- 
tically coincide with those observed in the ZFC mode. 
This shows that in our thin plate-like mesas suppres- 
sion of Josephson coupling by switching into the resistive 
state irreversibly produces pancake disorder comparable 
to that produced by ramping the magnetic field, i.e. vor- 
tices do not change their positions when one goes from 
the resistive state back into the superconducting state. 

In Fig. |j| we present our principal result showing a 
typical dependence of QP conductivity as a function of 
magnetic field. The conductivity rises steeply with mag- 
netic field, reaches a distinct plateau, and then rises more 
gradually with field at higher fields. Below we explain our 
experimental data using theoretical results |ll]] for a q (B) 
and an estimate for J C (B,T) in small mesas. 

Let us outline the main theoretical predictions |ll[] for 
QP conductivity <J q (B) in the mixed state. QP tunneling 
in a d-wave layered superconductor in the presence of vor- 
tices is determined by two competing mechanisms: a q is 
proportional to the QP DOS and inversely proportional 
to the effective scattering rate for interlayer tunneling, 
7 C . In-plane supercurrents around vortices lead to a DOS 
increase near gap nodes proportional to the averaged 
Doppler shift in the QP spectrum, E H ~ vf{B /^o) 1 ^ 2 
at 7 <C Eh "C Ao , where 7 is the impurity induced effec- 
tive in-plane scattering rate and Ao is the gap amplitude. 
In the framework of the specular tunneling model, when 
in-plane momentum is conserved at tunneling, a q will in- 
crease strongly with B due to DOS enhancement when 
vortices are correlated along the c-axis. If c-axis corre- 
lation of pancakes is absent, the effective scattering rate 
7 C increases with B due to different Doppler shifts at 
equivalent points in adjacent layers, suppressing tunnel- 
ing. There is a characteristic field, B 1 — $o7 2 /h 2 Vp, es- 
timated as B 1 ~ 0.2 — 0.6 T, which corresponds roughly 
to the field at which the variation of Doppler shift be- 
comes comparable to the scattering rate 7. The lat- 
ter may be estimated from the temperature dependence 
a q (0,T) = cr 9 (0, 0)(1 + cT 2 ), where c = ir 2 /l8~f 2 . Values 
of 7 for mesas studied are presented in the Table I. It 
was predicted that in the quasiclassical approach, valid 
at Eh <C Ao, when one accounts for Doppler shifts only, 
one gets 7 C ~ Eh in fields B > £? 7 . In this field range 
increase in scattering caused by Doppler shift variation 
in adjacent layers compensates in a q {H) the increase in 
DOS, providing a nearly field independent a q (H). Cor- 
rections to this approach lead to a weak linear increase 
of QP conductivity at B 1 <C B <C Bq, 



were calculated for a 2D vortex liquid depending on the 
effective temperature of vortex disorder T e / / . 



a q (B)/a q (0)=C 1 + B/B o , 



(1) 
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FIG. 4. For mesas [ml (□), m2 (o) and m3 (v)] the nor- 
malized o q (top) and (J c ) (bottom) are plotted vs. H. Three 
different regimes with respect to H are marked by verti- 
cal dotted lines. The theoretical curves of cr q (B) for cor- 
related (dashed) and c-axis-uncorrelated (dash-dotted) vor- 
tex systems are shown for B 1 = 0.6 T, Ao = 25 meV, and 
T ef} =0.06 ($ 2 s/8ir 2 \ 2 ab ) 0. 

The field dependence of the Josephson critical current 
in the vortex glass state for H || c comes from the depen- 
dence of J C {B, T) = J c (0, T)(cos (r)) on the gauge- 
invariant phase difference ip n ^ n+ i(r) induced by pancakes 
when they are misaligned in the layers n and n + 1 due to 
random pinning. Correlations of pancakes in neighboring 
layers, leading to lower energy, are induced by magnetic 
coupling of pancakes and by Josephson interlayer cou- 
pling; both tend to align vortices along the c-axis. Maxi- 
mum J c in the vortex glass phase may be achieved by ad- 
justment of pancake positions to minimize both Joseph- 
son and pinning energy [ [j"6[ . This effect is negligible for 
the ZFC mode and, as shown above, for the FC mode 
after switching from the resistive state where Josephson 
coupling is absent. For the critical current J c (B,0) in 
mesas with the size I < Aj/a, the relations 



(J?(B,0)) 



drdr' 



(exp[i^„, n+ i(r) - iip n>n+ i{r')]), (2) 



where Bq w 20 T. Here C\ depends on in-plane pan- 
cake ordering. The values of C\ between 1.07 and 1.22 



J 2 (0,0) J 

i 

hold |fl3|| . Here a = (^q/B) 1 / 2 is the intervortex distance, 
Xj is the Josephson length, ~ 1 [im in typical mesas, 
4> v (r) is the polar angle of the point r, and r ni is the co- 
ordinate of pancake i in the layer n. In Eq. (^) the inte- 
grand, (exp[iip nt n + i(r)-i<p nin+1 (r')]), drops with |r-r'| 
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on the scale a at a <C L. Then we estimate (J C (B,Q)} < 
{Jl(B,$)) x / 2 w J c (0,0)($o/-BL 2 ) 1/2 . Thermal fluctua- 
tions cause jumps into the resistive state at J < J C (B, 0) 
and diminish additionally the average critical current. 
When the Josephson energy, & (J C (B, 0)}L 2 /2ttc, be- 
comes less than » 10T (see Eq. (6.18) in Ref. Q), ther- 
mal fluctuations practically eliminate the superconduct- 
ing part of I-V curves. For our mesas this occurs in fields 
> 0.01 T at T = 4 K. 

Now we explain our experimental results using these 
theoretical results. In Fig. [| we plot the low-field part of 
the <J q {H) dependence together with extracted J c ,max{H) 
for three different samples. Shown here are also theo- 
retical curves calculated for cr q (B) in the case of c-axis- 
correlated and uncorrelated vortex states H] , indicating 
clearly that c-axis disorder is responsible for strong sup- 
pression of <j q at high fields. The experimentally observed 
<J q {H) can be divided into three segments as magnetic 
field increases as shown in Fig. 4. These segments cor- 
respond to appropriate parts of the field dependence of 
the critical current density, J c ,max, obtained in the ZFC 
mode, presented in the lower part of Fig. 4. At very 
low fields, H < 5 mT, (Segment I) there is practically 
no change in a q and in J c .max- This regime corresponds 
to the Meissner state of the sample. Note that, despite 
the suppression of supercurrent along the c-axis by ap- 
plied transport current in the resistive state, the in-plane 
Meissner currents persist, preventing pancake vortex en- 
trance into the sample. Only when the magnetic field 
H exceeds ~ 5 mT do pancakes enter into the sample. 
In the intermediate field range, 5 < H < 40 mT, (seg- 
ment II) the experimental results for a q (H) lay slightly 
above the curve for uncorrelated vortices but well below 
the curve for the c-axis-correlated vortex state. This in- 
dicates that weak c-axis vortex correlations are present 
in this field range. Weakly correlated pancakes affect J c 
much more strongly than a q , because increase in B leads 
to decrease of J c ( the fluctuations of the phase difference 
increase and (cosy„. n +i(r)) drops), while for a q the in- 
crease of DOS is compensated partly by an increase of 
7 C . As the field continues to increase above 40 mT (seg- 
ment III) increase of a q with H is slowing down due to a 
stronger decrease of c-axis correlations and, respectively, 
stronger increase of 7 C . This leads to formation of the 
plateau (or even dip) in the range 40 < H < 300 mT. 

This scenario is confirmed by the behavior of 
J c ,max(H). A sharp drop in c-axis correlations above 
40 mT was observed in Bi-2212 single crystals as the sec- 
ond peak effect. It was shown by neutron scattering |17| ] 
and JPR frequency study |n| that here the vortex sys- 
tem decomposes into a c-axis-uncorrelated vortex glass. 
Accordingly, the c-axis critical current density, J cm ax, 
drops rapidly, as seen clearly in the insert in Fig. [l]. At 
fields above 0.5 T the dependence <J q (H) changes again, 
switching to the linear growth, observed and discussed 
in our previous work [pfl. We obtain B ~ 20 T and 



C\ ~ 1.04. We see that C\ is smaller than that calcu- 
lated for the vortex liquid, and this means that in-plane 
ordering in the vortex glass is better than in the liquid. 

To conclude, we observed a plateau in a q {H), which 
formed as a result of the compensation of the increase 
of QP DOS due to Volovik effect by the increase of the 
effective scattering rate for tunneling due to c-axis pan- 
cake disorder. The field dependence of the QP interlayer 
conductivity is sensitive to the structure of vortex state. 
We found that in the vortex glass adjustment of vortex 
positions to the Josephson interlayer coupling depends on 
the way this vortex state was achieved. Critical current 
in the FC mode is ~ 2 times larger than that obtained in 
the ZFC mode or in the FC mode after switching from 
the resistive state. The c£-wave pairing model in the clean 
limit with resonant intralayer scattering and significant 
contribution of specular interlayer tunneling provides an 
explanation for the c-axis electrical transport in the vor- 
tex state of Bi-2212 crystals at low temperatures and low 
magnetic fields. 
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